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Abstract. We present X-ray observations of the interacting 
pair of galaxies NGC 4410a/b with the ROSAT HRI and PS PC. 
The ROSAT HRI images reveal a point-like source correspond- 
ing to NGC 4410a and an X-ray halo, extending 10" from the 
nucleus toward the southeast. The halo emission accounts for 
-1/3 (L x = 1.3xl0 41 erg s" 1 in the 0.1 - 2.4 keV ROSAT 
band) of the total X-ray emission detected from NGC 4410a. 
The spectrum of the total X-ray emission from NGC 4410a 
can be fitted at best with a two-component emission model, 
combining a Raymond-Smith spectrum with T = 10 7 K and a 
power-law spectrum (T = 2.2). The fraction of the thermal com- 
ponent to the total flux within this model is 35%, supporting the 
results of the HRI observation. The total unresolved X-ray lu- 
minosity detected with the PSPC amounts to 4x10 erg s~ 4 
in the ROSAT PSPC band (for D = 97 Mpc). 

A preferable explanation might be the interaction between 
the two galaxies causing a circumnuclear starburst around the 
active nucleus in the peculiar late type Sab galaxy NGC 4410a. 
As a cumulative effect of exploding stars a superbubble forms 
a bipolar outflow from the galactic disk with an expansion time 
of —8 Myr. The central source injects mechanical energy at a 
constant rate of a few times 10 42 erg s~ 4 into the superbubble. 

A HST WFPC2 image decovers that NGC 4410a is seen 
almost face-on so that only the approaching outflow is visible. 
The HRI contours reveal an elongation of this outflow indica- 
tive for a faint X-ray ridge below the detection limit (Lx < 
1 .4 x 1 39 erg s ~ 1 ) toward the neighbouring galaxy NGC 44 1 0b 
caused by its tidal interaction. NGC 4410b houses only a faint 
X-ray source (Lx < 3.8 x 10 39 erg s~ 4 ). 

Key words: Galaxies: interactions - Galaxies: active - Galax- 
ies: starburst - Galaxies: individual: NGC 4410 - X-rays: 
galaxies 



1. Introduction 

The pair of interacting galaxies NGC 4410a/b belongs to a 
group of 1 1 members (Hummel et al. 1986, hereafter HKG86) 
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which are located behind the Virgo cluster. It consists of a pe- 
culiar Sab (NGC 4410a) (Thuan & Sauvage 1992) and an E 
galaxy (NGC 4410b) (HKG86) located in east- west direction 
and separated by 18'.' 7 (8.8 kpc at a distance of 97 Mpc) (Maz- 
zarella & Boroson 1993, hereafter MB93). NGC 4410a is lo- 
cated at RA = 12 h 26 m 27!9, Dec = +09° 01' 18" (J2000). 
Spectral analysis of the nuclei of both components led to the 
classification as two LINERs (Mazzarella et al. 1991; Bicay et 
al. 1995; Thuan & Sauvage 1992). In the eastern component 
NGC 4410b a supernova type I has been detected in 1965 (SN 
1965 A) from which Turatto et al. (1989) derived a distance 
of 139 Mpc. The system has a radial velocity of about 7300 
km s" 1 (Mazzarella et al. 1991; Batuski et al. 1992; Thuan & 
Sauvage 1992; MB93). With a Hubble constant of 75 km s" 1 
Mpc -1 this leads to a distance of 97 Mpc. In this paper we ap- 
ply the latter value to all distance-dependent parameters, like 
e.g. luminosities etc. A distance of 97 Mpc for NGC 4410 re- 
sults in an absolute length scale of 470 pc arcsec -1 . 

HKG86 found a spatial coincidence of the western optical 
component NGC 4410a with a radio point source of luminosity 
I/r psIO 39 erg s _1 embedded in a strong, extended radio source 
around NGC 4410a with a total radio luminosity L R «7x 10 40 
erg s _1 . 

The phenomenon of interaction between galaxies is closely 
related to the occurrence of starbursts (SBs): A large fraction 
of interacting galaxies (~70%) exhibits typical characteristica 
of SBs (Bushouse 1986) and vice versa. It is also striking that 
the majority of infrared-bright galaxies shows evidence for re- 
cent interactions as indicated by the presence of close neigh- 
bours, or by their disturbed morphology and tidal tails (Joseph 
et al. 1984; Lonsdale et al. 1984; Telesco 1988). The frac- 
tion of mergers increases drastically in the infrared luminosity 
range L IR from 10 10 L Q (12%) to 1O 12 L (95%) (Sanders & 
Mirabel 1996). Although a large IR luminosity is not necessar- 
ily a tracer for enhanced star formation, it is one of the typical 
features for SB galaxies. 

Norman & Scoville (1988) demonstrated that gas in- 
flow from the galactic disk toward the central region, as re- 
quired for building a central SB, will be triggered by a non- 
axialsymmetric gravitational potential like in interacting galax- 
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ies. They argued that perturbed orbits of gas clouds caused 
by an encounter lead to enhanced cloud collisions within the 
galactic central region and trigger the formation of massive 
stars. An alternative model for triggering a central SB was 
given by Jog & Das (1992) and Jog & Solomon (1992). The 
infall of giant molecular clouds (GMCs) into an intercloud 
medium with a higher mean pressure in the central region 
drives a radiative shock into the GMCs and ignites the SB. 

The inflow of gas from the galactic disk is not only required 
for a central SB but also for fuelling an AGN. It is still under de- 
bate whether SBs are progenitors of AGN (Norman & Scoville 
1988) or whether both are different physical processes. Weed- 
man (1983) proposed that if a large number of massive stars 
form fast in a small central volume, the compact stellar rem- 
nants from these could act as accretors. Also the stellar dy- 
namical merging of a dense cluster of massive stellar remnants 
would plausibly form a blackhole nucleus (Rees 1984). 

The nuclei of several nearby galaxies, like NGC 1068, 
NGC 1097, and NGC 7469, can be resolved into a central AGN 
and a circumnuclear SB ring (Keel 1985; Perez-Olea & Colina 
1996, hereafter PC96). Some host galaxies of quasars also ex- 
hibit evidence for SBs, such as 3C48 (Stockton 1990). 

Assuming an evolution from starbursts to AGNs during the 
interaction of galaxies one would expect more AGNs with the 
age of the merger process. Recent observations of ultralumi- 
nous IR galaxies with the infrared satellite ISO however do 
not show any obvious tendency of an increasing fraction of 
AGNs within interacting galaxies with advanced merging pro- 
cess (Lutz et al. 1998; Genzel et al. 1998). Half of the observed 
galaxies reveal both an AGN and starburst activity. It seems 
more likely that more local and shorter term conditions like 
time-dependent compression of the circumnuclear interstellar 
gas, the accretion rate onto the central black hole, and the as- 
sociate radiation efficiency determine AGN or starburst domi- 
nated luminosities. 

In this paper we present the spectral and imaging X-ray 
properties of NGC 4410 observed by the ROSAT Position Sen- 
sitive Proportional Counter (PSPC) and High Resolution Im- 
ager (HRI), respectively. We are able to resolve the two optical 
components in the ROSAT HRI image. The paper is structured 
as follows: In the next section we are considering the X-ray ob- 
servations in general, detecting and trying to identify the emer- 
gent X-ray sources in the field and discussing the data reduc- 
tion. In Sect. 3 we concentrate on NGC 4410 presenting the 
HRI results at first before we study the spectral flux distribu- 
tion of the PSPC data with particular concern of comparison 
with a sequence of X-ray radiation models. The results are then 
discussed in Sect. 4 making use of a HST WFPC2 image of 
NGC 4410a for a geometrical consideration of its inclination. 

2. Observations and data reduction 

The data presented in this paper were taken with the HRI and 
the PSPC detectors on board of the X-ray satellite ROSAT. This 
X-ray telescope is operating in the energy range of 0. 1 to 2.4 
keV. The spatial resolutions of the HRI is l'/7. The point spread 
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Fig. 1. Central 42' of the PSPC field of view centered around 
NGC 4410 (bright central source) with an exposure time of 
23.4 ksec. The image is smoothed with a Gaussian filter (20" 
FWHM). The scale of the image is slightly different from the 
one in Fig. 2. North is up, east is to the left. 



function (PSF) of this detector at the optical axis in combi- 
nation with the telescope is 3". The PSPC has a PSF of 25". 
The two detectors have a field of view of 38' and 2°, respec- 
tively. For details concerning ROSAT and its instruments see 
the ROSAT User's Handbook (Briel et al. 1996). 

NGC 4410 was observed on June 28-30, 1993 with the 
ROSAT PSPC detector for a total effective exposure time of 
23.4 ksec. The total number of background-subtracted counts 
from the central source associated with NGC 4410 is 870±31. 
Figure 1 shows the central 42' of the PSPC field of view. For 
the spectral analysis of the source we used the software pack- 
age IDL. The source photons were extracted from a circular 
area of 92" around the central source and corrected for tele- 
scope vignetting and detector dead-time. The spectrum of the 
source was binned according to a signal-to-noise ratio of 8. 
For the background correction we selected three uncontam- 
inated circular areas close to NGC 4410 with radii of 86", 
101" and 141". The background contributes an X-ray flux of 
(1.36±0.27)xl0~ 2 cts arcsec~ 2 in the ROSAT bandpass. In 
order to fit a model spectrum to the PSPC data we used the X- 
ray spectral-fitting software package XSPEC (Arnaud 1996). 

Our ROSAT HRI observations have been taken between 
June 27 and July 1, 1995 with a total effective exposure time 
of 36.7 ksec which led to 152±24 counts for the central source 
associated with NGC 4410. The background X-ray flux amouts 
to (1.6±l.l)xl0~ 2 cts arcsec~ 2 . 
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Fig. 2. ROSAT HRI X-ray sources around NGC 4410 (source 
no. 8) detected by the maximum likelihood method and with a 
threshold of 5a (a = 3. Ox 10~ 7 cts s _1 arcsec~ 2 ), plotted over 
an optical image taken from the Digitized Palomar Observatory 
Sky Survey. The size of the field corresponds to the HRI field 
of view. (Note the different scale compared to Fig. 1.) North is 
up, east is to the left. 

2.1. Source detection 

Ten X-ray sources above a limit of 5 a are found within the 
field of view of the HRI detector using a standard source- 
detection algorithm in the EXSAS X-ray analysis software 
package (Zimmermann et al. 1997) (Fig. 2). In Table 1 the de- 
tected sources are listed with position, background subtracted 
count rates and optical identifications. 

From these detected X-ray sources in the HRI image, four 
coincidences with optical objects can be found (see Table 1 
and Fig. 2). The X-ray source no. 8 can be clearly identified 
with NGC 4410a, while an extended X-ray source (no. 10, 44" 
FWHM) belongs to the galaxy cluster Abell 1541. In addition, 
the source no. 1 coincides with the QSO 1224+095 and no. 4 
is still unidentified in the optical. In order to investigate fur- 
ther identifications of X-ray sources with optical objects we de- 
crease the detection limit to 3a. This leads to 32 sources in the 
HRI image with 10 counterparts in the optical. Three of them 
are identified as the Abell galaxy cluster 1541 (extended), the 
QSO 1224+095 (1" offset) and the K7-star G60-2 (8" offset). 
The offsets between optical and X-ray sources do not indicate 
any systematical shift or rotation, requiring a position correc- 
tion of the superposition of the two images at this point. 

3. Results 



3.1. X-ray imaging 

Using the coordinates of the ROSAT HRI pointing the overlay 
of X-ray and optical image reveals a shift between the optical 
(NGC 4410a) and X-ray maximum of 4". Strikingly, an addi- 
tional X-ray source, that is visible at the 3a contour level at a 
distance of ~48" from the central X-ray maximum to the east, 
has the same displacement from the maximum of the eastern 
optical component NGC 4410b. It contains about 1% of the 
count rate of the central source in the HRI image. From the 
detection limit of 3 a = 9.0 x 10~ 7 cts s _1 arcsec~ 2 in the HRI 
image and the energy conversion factor (ECF) of 6.7 x 10 9 cts 
cm 2 erg -1 for the HRI and a power-law spectrum with T = 2.4 
from the ROSAT User's Handbook (Briel et al. 1996) we get 
an upper limit of 3.8 x 10 39 erg s _1 for the X-ray luminosity of 
NGC 4410b. 

Furthermore, radio observations of NGC 4410 present a 
similar emission feature as the X-ray contours. The maximum 
radio contours correspond spatially to the optical maximum of 
NGC 4410a and the lower emission levels are elongated to- 
ward the southeast (HKG86). To check whether the eastern 
X-ray maximum has to coincide with the optical structure of 
NGC 4410, we compare its X-ray luminosity with the blue lu- 
minosity L B of NGC 4410b. With m B = 15.28 (MB93) we de- 
rive a blue luminosity L B = 1.2 xlO 10 L Q for NGC 4410b. 
Brown & Bregman (1998) analysed a sample of 34 early-type 
galaxies observed with ROSAT PSPC and HRI and studied the 
X-ray-to-optical distribution. The estimated X-ray luminosity 
of ~ 10 39 erg s -1 for NGC 4410b fits well into the observed 
distribution of the sample. This obtrudes to superpose the X-ray 
maxima of both components NGC 4410a and b to the optical 
image of the pair. This fine tuning requirement is not detectable 
and therefore not necessary on the scale of the full HRI image 
(Fig. 2) as pointed out in Sect. 2.1. 

Figure 3 shows the position corrected contour plot of the 
central 2' of the HRI image with the X-ray source overlaid onto 
an optical image of NGC 4410 of the Digitized Palomar sky 
survey plates. The HRI image was smoothed with a Gaussian 
filter FWHM of 3". 

The contours of the HRI image indicate a radial extension 
from the maximum of NGC 4410a toward the southeast by 3" 
to 10". Since this elongation is close to the resolution limit it 
has to be checked whether it is real or artificial. Morse (1994) 
found out that in some ROSAT HRI observations the contours 
of point-like X-ray sources are elongated over a scale of ^5" 
to 10" because of errors in the attitude correction as the satel- 
lite was wobbled during the observation. On the other hand, 
it is striking that the extended radio emission of NGC 4410a 
(HKG86) agrees in its direction with the X-ray contours. More- 
over, the elongation of the X-ray contours is not symmetrical 
with respect to the elongation axis, as one would expect if it 
is caused by the wobbling satellite, similar to the point-like 
sources in the paper by Morse (1994). Another possible error 
namely that the elongation is caused by an incorrect superposi- 
tion of the different used observation intervalls (OBIs) was also 
checked by us but didn't explain the effect. 



4 D. Tschoke et al.: ROSAT X-ray observations of NGC 4410 



Table 1. ROSAT HRI detections in the NGC 4410 frame with a maximum likelihood threshold of 5a. 
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Fig. 3. Overlay of the HRI X-ray contours onto the optical im- 
age of NGC 4410 taken from the Digitized Palomar sky survey. 
The X-ray image was smoothed with a Gaussian filter FWHM 
of 3". The contour levels are at 3, 5, 7, 10, 15, 20, 30 and 
35ct above the mean background with a = 3.0xl0~ 7 cts s -1 
arcsec -2 . The relative position of optical image and X-ray con- 
tours was corrected as explained in Sect. 3.1. North is up, east 
is to the left. 
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Fig. 4. Residual X-ray flux after subtracting a point-like source 
with Gaussian distribution (FWHM 3.5") at the position of the 
maximum of the X-ray flux from the original total flux image 
(see Fig. 5). The cross marks the position of the maximum 
and coincides with the nucleus of NGC 4410a. The contours 
are 3, 5, 7, 10 and 13cr above the mean background with a = 
3. Ox 10~ 7 cts s _1 arcsec -2 . North is up, east is to the left. 



We fit a Gaussian distribution with FWHM 3'.'5 (spatial res- 
olution of HRI = 1'.'7, image smoothed with Gaussian filter 
FWHM 3" corresponding to the HRI PSF) to the point-like 
source and subtract it from the original image. The result is 
plotted in the Figs. 4 and 5. The residual reveals a ~10"x5" 
emission feature 7" to the southeast of the X-ray maximum vis- 



ible up to 13cr above the background. As expected, the wings 
of the Gaussian distribution in the Figs. 5b-f fit very well to 
the X-ray source except in the direction of deformation of the 
X-ray contours. 

As a working hypothesis which will be supported by the 
finding of a thermal emission (see Sect. 3.2) we assume that 
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Fig. 5. a) Positions of the five sections (dashed lines) through the HRI image of NGC 4410. The contours are the same as in 
Fig. 3. North is up, east to the left. The following five plots show the corresponding sections through the X-ray maximum of 
the HRI image. The dashed lines in all five section plots represent the fitted Gaussian distribution with a FWHM 3'.'5 to the 
point-like source of NGC 4410a. b) Section 1, east-to-west. The local maximum 18" east of the X-ray peak spatially corresponds 
to NGC 4410b. c) Section 2, southeast-to-northwest, d) Section 3, southeast-to-northwest, e) Section 4, southwest-to-northeast, 
f) Section 5, southwest-to-northeast. In each plot the Gaussian distributed compact source (dashed lines) fits well to the HRI 
contours except in the direction of deformation of the X-ray contours. 




Fig. 6. Best single-component model fits to the observed (0.1-2.4) keV X-ray spectrum of NGC 4410, PO (left) and BS (right). 
The channel counts are binned to a signal-to-noise ratio of 8. The parameter values for these models are shown in Table 2. The 
lower graphs show the residuals of the fits. 



the extended X-ray contours represent a superbubble originat- 
ing from type II SNe in the starburst region, and expanding 
into the halo. Since its elongation is too small to resemble the 
edge-on view of an outflow we further assume that the residual 
emission in Fig. 5 represents the front part of a bipolar out- 
flow from the central region of the galactic disk. To verify this 
assumption, we estimate the ratio of the integrated flux from 
the residual image with the total flux of the original image. For 
that we integrate the flux within a radius of 13" from the maxi- 
mum in both the residual image and the original one. We derive 
a ratio of 0.34±0.05 (52±8 counts). The error is estimated by 
varying the integration radius around the source by 2" and from 
the mean fluctuation of the background flux. This result is sup- 
ported by the spectral analysis of PSPC data where a flux of Lx 
= 1 .3 x 10 41 erg s _1 for the thermal emission component is de- 
rived as part of the combined thermal and power-law spectrum 
(see Sect. 3.2). However, since the fraction of the superbubble 
within the point-like source is neglected, this estimate is only a 
lower limit. 

An X-ray ridge toward the optical position of NGC4410b 
is discernable at the 3a level, which corresponds to an up- 
per count level of 3. Ox 10 -7 cts s^ 1 arcsec~ 2 . Using the HRI 
ECF of 1.8x 10 12 cts cm 2 erg -1 for a 1 keV Raymond-Smith 



plasma with log iVn = 20.0 from the ROSAT User's Handbook 
(Briel et al. 1996) one gets an upper limit for the X-ray lumi- 
nosity of 1.4 xlO 39 erg s _1 . Two regions of enhanced X-ray 
emission are visible 15" and 22" south of the maximum X-ray 
source that were not contained in our source list with a 3a con- 
fidence level. Depending on the spectral model and therefore 
on the ECF the upper limit for the X-ray flux of these sources 
is2.25xlO- 16 /ECFergs- 1 cm" 2 . 

3.2. Spectral analysis 

In addition to the spatial information of emission components 
we got from the HRI data, we now want to decompose these 
components by using the spectral information obtained from 
the PSPC data. With the PSF of 25" of the PSPC we get no 
satisfactory spatial resolution. 

The 870 background corrected counts of the central source 
detected with the PSPC in the energy range of 0. 1 -2.4 ke V were 
binned in order to achieve a constant signal-to-noise ratio of 
8. To fit the observed PSPC X-ray spectrum we apply differ- 
ent models of absorbed emission spectra including a thermal 
Bremsstrahlung spectrum (BS), the spectrum of an optically 
thin thermal plasma in collisional equilibrium (RS) as calcu- 
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Table 2. Spectral fits to the PSPC data of NGC 4410 and the derived model parameters. For the determination of the X-ray 
luminosity a distance of 97 Mpc is assumed. 
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- 


2.6 
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76 
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0.7 
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fixed to Galactic foreground value 
separate absorption for each component 



Col.(l) — Emission models abbreviated as: BS = thermal Bremsstrahlung, RS = Raymond-Smith, PO = power-law. The specification 
in brackets indicates the used element abundances for the hot thermal plasma, solar (sol) or abundances of a supernova calculated by Woosley 
and Weaver (ww) and integrated over all stars between 10 and 100 Mq for a Salpeter IMF. 
Col. (2) — column density in units of 10 20 cm~ 2 . 
Col. (3) — plasma temperature in units of 10 6 K. 
Col. (4) — photon index. 

Col. (5) — scaling factor: for BS in units of (10~ 19 /(47rD 2 ))J n c nidV, n c ,n-a = electron and ion densities (cm -3 ); for RS in units of 

(10~ 20 /(4ttL> 2 )) / n c n H dV\ n e , n H = electron and H densities (cm -3 ); for PO in units of 10~ 5 photons keV -1 cnT 2 s _1 at 1 keV. 

Col. (6) — reduced \ 2 - 

Col. (7) — degrees of freedom. 

Col. (8) — X-ray flux in units of 10~ 13 erg cm~ 2 s _1 , corrected for the fitted absorption. Values in brackets give the contribution of the thermal 
component. 

Col. (9) — X-ray luminosity in units of 10 41 erg s^ 1 . Values in brackets give the contribution of the thermal component. 



lated by Raymond & Smith (1977) and a power-law spectrum 
(PO) with a photon energy distribution oc E~ r , as well as sev- 
eral combinations of those. The absorption of the spectra de- 
pends on the hydrogen column density iVn between the source 
and the observer. As a first-order estimate and as a lower limit 
of iVn we use the Galactic foreground value in the direction 
of NGC 4410. This value for N H amounts to 1.73 x 10 20 cm" 2 
(Dickey & Lockman 1990). 



The model fits for the X-ray emission to the observation 
are shown in Table 2. Usually the RS model is applied under 
the assumption of solar abundances (Anders & Grevesse 1989). 
Reasonably, this is not an ideal approach because expanding 
SN ejecta should contain the nucleosynthesis products of their 
massive stellar progenitors. For this reason we use two different 
sets of chemical abundances. Beside the classical model with 
solar abundances we fit one with enhanced abundances of SNe 
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99%, 90% and 68% confidence. 
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according to the yields by Woosley & Weaver (1985), but inte- 
grated over all massive stars between 10 and 100 Mq and with 
a Salpeter initial mass function (IMF). 

The main uncertainty to fit the data comes in by the value 
of the hydrogen column density Ah- In each model, except the 
combinations of RS+PO and the single component models BS 
and PO, the free parameter A^h turns out to become much too 
low for the best fit. For the single thermal model RS it is not 
possible to fit both, the soft and the hard part of the spectrum 
reasonably without an unrealistically low value of A^h- 

This is also true for the combined two temperature plasma 
model RS+RS(ww). The best fit results in an unrealistically 
low Ayn value of 0.35 xlO 20 cm -2 . For the same model 



> 

CD 



\ 0.010 - 



c 
=1 
o 
u 



0.001 

0.6 
0.4 

0.2 

0.0 

-0.2 
-0.4 




0.5 



1.0 1.5 2.0 

Energy/kev 



2.5 



Fig. 9. Observed X-ray spectrum of NGC 4410 and best fit 
(solid line) with a combined RS and PO model(RS+PO(ww)). 
The parameter values for this model are shown in Table 2. 
The long-dashed line shows the spectrum of the thermal com- 
ponent. The short-dashed line shows the contribution of the 
power-law. The contribution of the RS component to the to- 
tal flux is 35%. The residuals of this fit are plotted in the lower 
box. 



with solar abundances RS+RS(sol), the best fit yields Ayn = 
0.45 x 10 20 cm" 2 . In addition, the plasma temperature of more 
than 3xl0 7 K for the second component is abnormally high 
for halo gas. This temperature may only be reached for gas in 
the central regions of starburst galaxies, for example NGC 253 
(Ptak et al. 1997) and M82 (Mathews & Doane 1994). In the 
case of fixing the column density to the Galactic value, the re- 
sulting plasma temperatures of 1.7x 10 6 K and 8.8 x 10 6 K for 
RS+RS(sol) and 0.9xl0 6 K and 9.3xl0 6 K for RS+RS(ww) 
are much more reasonable. But both fits to the observed spec- 
tra are much worse (see Table 2). 

In the case of fitting the observation with a single PO the 
column density of 3.4 x 10 20 cm~ 2 becomes plausible because 
of an expected additional absorption within NGC 4410. The 
photon index of 2.39±0.09 is in agreement with the results 
from Turner et al. (1993) and Mulchaey et al. (1993), who 
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model (Fig. 9). The levels represent 99%, 90% and 68% confi- 
dence. 



found a mean value for Seyfert 1 and Seyfert 2 galaxies of V = 
2.4. Comparing all these single-component models the model 
PO and BS provide the best fits to the observation. The fits are 
plotted in Fig. 6. 

The best fit by far for the whole observed spectrum can be 
reached by applying a two-component model which combines 
a thermal and a power-law spectrum (RS+PO) with a photon in- 
dex of 2.17±0.10 and a plasma temperature of (9.3±1.0)x 10 6 
K (see Fig. 9). Nevertheless, the fit becomes worse for apply- 
ing a higher value for Nh than the Galactic one. The combined 
PO+RS fit yields the best result with an absorption corrected 
luminosity of (3.91±0.55)x 10 41 erg s _1 for the assumed dis- 
tance of 97 Mpc. 

Under the assumption, that the nonthermal emission orig- 
inates from a central active nucleus heavily obscured by dust 
and gas, and that hot plasma has expanded into the outer parts 
of NGC 4410a, we fit the RS+PO model with different column 
density values for each component. One would expect a column 
density -/Vh,rs ~ -ZVn.Gai for the less obscured plasma outflow 
and a much higher -/Vh,po f° r the AGN because of the intrin- 
sic absorption within the NGC 4410 nucleus. While we apply 
A\h,rs = A\H,Gai = 1.73 xlO 20 cm" 2 for the RS component, 
TVh.po for the PO component is set to 3.1 x 10 21 cm~ 2 . This 
yields the best fit for a lower plasma temperature of 1 .0 x 10 6 K 
and a very steep power-law for the central AGN with T = 4.24 
(see Fig. 1 1). In the very soft range (0.1-0.3 keV) the spectrum 
is determined by the intrisically unobscured RS component. 
Raising the absorption for the PO component prohibits its con- 
tribution to the soft spectral range, which means that the total 
X-ray flux there has to originate from the hot ionized gas. As a 
consequence the temperature of the plasma is almost one order 
of magnitude lower than in the RS+PO model with only one 
common absorption. The hard spectrum >0.5 keV originates 
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Fig. 11. Fit of the RS+PO(ww) model with individual column 
densities for each component and residuals of the fit. The dif- 
ferent lines represent the same components as in Fig. 9. As 
expected the soft part of the spectrum is dominated by the un- 
absorbed thermal emission from the hot plasma. The hard spec- 
trum >0.5 keV originates from the power-law component. Yet 
the spectrum can only be fitted with much too steep photon 
index of more than 4. 

completely from the nonthermal component. Nevertheless the 
photon index is unplausibly steep. 

4. Discussion 

4.1. The X-ray halo 

The HRI image strongly suggests a twofold spatial correlation 
between the X-ray and the radio emission: we have a compact 
central emission peak apparently coinciding with NGC 4410a 
and a diffuse halo slightly elongated about 10" toward the 
southeast. Subtraction of the central point-like source from the 
total X-ray emission enables us to unveil a clear extention in 
this direction. Because of the asymmetrical elongation and the 
spatial correlation in the X-ray and radio flux we are convinced 
not to be mistaken by an artefact. Furthermore, the flux ratio 
of 0.34±0.05 from this feature compared to the total flux cor- 
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responds to the best fitting model for the observed spectrum, a 
two-component model, combining a thermal Raymond-Smith 
spectrum and a power-law spectrum. As a hypothesis we as- 
sume that the thermal emission originates from an expanding 
X-ray halo gas with a plasma temperature of a few 10 6 K whose 
expulsion is driven by a central engine (SB + AGN). 

The SE elongation of the emission contours could be 
caused by two effects: 

(1) If the pec Sab galaxy NGC 4410a is seen at moderate in- 
clination, a circumnuclear superbubble or a bipolar outflow is 
visible in soft X-ray which expands vertically to both sides out 
of the galactic disk only the part in front of the Hi disk. X-ray 
emission in the ROSAT band originating from gas outflowing 
into the reverse halo hemisphere is absorbed by cold gas in the 
disk of NGC 4410. This fact also accounts for the low N H that 
is required from the best model fit and is lacking of any addi- 
tional contribution to the Galactic foreground value. 

(2) The X-ray halo gas is exposed to the tidal force of the com- 
panion galaxy NGC 4410b. 

In order to decide which of the two effects causes the de- 
formation of the X-ray halo, information about the inclination 
angle of NGC 4410a is necessary. Fig. 12 shows the Hubble 
Space Telescope (HST) WFPC2 image of the central part of 
NGC 4410a in the V (5843 A) broad band filter taken from the 
Hubble Space Telescope archive. One can discern the nucleus, 
a bright emission region at the lower left, and a line of addi- 
tional smaller bright knots at the left. At the right-hand side of 
the nucleus a long tidal arm probably caused by the interaction 
is clearly visible. Since the morphology in this image leads to 
a very low inclination, we prefer the second explanation of an 
interaction-induced deformation of the X-ray halo. However, 
the HST image suggests an additional possibility, namely, that 
the bright emission region at the lower left is a starburst region, 
and that the extended part of the X-ray contours originates from 
this region. 

The X-ray contours reveal a faint ridge toward NGC 4410b 
(L>x < 1.4xl0 39 erg s" 1 ) which has only a very faint corre- 
sponding X-ray counterpart centrally of NGC 4410b (Lx < 
3.8 x 10 39 erg s~ 4 ) as indication of any activity in or around the 
nucleus of NGC 4410b. The non-detection of a similar radio 
emission admits at best a radioquiet AGN. 

4.2. Physical conditions of the X-ray emission 

As one can see in Table 2 the derived values for the X-ray lu- 
minosity are all of the same order of 3^1-xlO 41 erg s _1 and 
independent of the applied model. From this point of view it is 
not possible to distinguish between the different models. The 
two most acceptable models with reasonable values for Nh 
are the two-component combination RS+PO and the single- 
component PO model, respectively. 

Although the photon index of T = 2.39±0.09 of the sin- 
gle PO model is in good agreement with observed values for 
low luminousity AGNs (Turner et al. 1993; Mulchaey 1993), 
it is obvious that the observed X-ray spectrum of NGC 4410 
cannot be fitted properly by a single PO model but requires an 



additional component as can be seen in Table 2 and by com- 
paring Fig. 6 with Fig. 9. The RS+PO model for the X-ray 
source leads to a luminosity of the thermal component Lx = 
(1.35±0.19)xl0 41 erg s" 1 (35% of the total X-ray luminos- 
ity). Nevertheless, the X-ray luminosity can vary by a factor of 
4 because of the uncertainty of the distance to the source. To 
compare the derived properties of NGC 4410 with theoretical 
calculations, we apply the following model: We assume a cen- 
tral SB to drive an expanding superbubble. The central source 
injects mechanical energy since a time t cxp at a constant rate 
L mec h by sequential supernova explosions into the bubble. The 
supersonic expansion leads to a shock front and is heating the 
ambient gas. Assuming a homogeneous density of this ambient 
gas would lead to a spherically expanding bubble, as calculated 
by MacLow & McCray (1988). 

This model, however, does not correspond to the realistic 
conditions of a SB in a disk galaxy. Because of the vertically 
decreasing gas density the superbubble forms a "chimney" and 
blows out of the galactic disk into the halo by a bipolar outflow 
(Norman & Ikeuchi 1989). The dynamics, X-ray emission and 
disk-halo interaction of these outflows are modelled under dif- 
ferent conditions of ambient disk gas density, halo gas density, 
energetic input into the superbubble, and expansion time of the 
superbubble in a number of papers (Tomisaka & Ikeuchi 1986; 
Tomisaka & Ikeuchi 1988; MacLow et al. 1989; Suchkov et al. 
1994, hereafter S94; Michaelis et al. 1996). S94 found that the 
soft X-ray emission of the suberbubble in the range of (0.1- 
2.2) keV primarily originates from the shock-heated halo and 
disk gas with temperatures of 10 6 - 10 7 K, rather than from 
the supernova material itself. Depending on the density of the 
disk and halo gas, the scale-height of the disk, and the energy 
deposition rate from the SB, the morphology of the bipolar out- 
flow varies in structural parameters: opening angle within the 
disk, radius of the chimney, vertical extension of the superbub- 
ble, and ratio of the two last properties. Independent of these 
different conditions, the models always achieve a bipolar mor- 
phology. S94 derived an analytic expression for the shell X-ray 
luminosity in the (0.1-2.2) keV band: 

L X = 9.7 xlO 40 (i 42 ) 3/5 (no.oi) 7/5 (t7) 9/5 erg s" 1 , (1) 

where L 42 is the mechanical input energy in units of 10 42 erg 
s~ 4 , no. oi is the particle density of the halo gas in units of 0.01 
cm~ 3 , and is the expansion time of the superbubble in units 
of 10 7 yr. 

From the HRI observation we can estimate a maximum ex- 
pansion of the outflowing gas into the halo of 10", scaling 
to 4.7 kpc at a distance of 97 Mpc. The models of Tomisaka 
& Ikeuchi (1986), Tomisaka & Ikeuchi (1988), MacLow et al. 
(1989) and S94 reveal chimney radii within the galactic disk 
between 200 pc and 600 pc. The much more extended X-ray 
emission observed in NGC 4410a can therefore only originate 
from the outflow into the less dense galactic halo where it can 
expand in each direction. This involves that the SB has to be old 
enough in order to account for the escape of the expanding su- 
perbubble from the disk into the halo. That means a lower limit 
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Fig. 12. Broad band image in the V (5843 A) filter of the nucleus of NGC 4410a with a logarithmic grayscale taken from the 
Hubble Space Telescope archive. As a consequence of the interaction, a chain of knots and dust, a bright emission region near 
the nucleus and a long tidal arm are visible. The image has a size of 34"x34". North is up, east is to the left. 



of t cxp > 4 Myr for the expansion time (Tomisaka & Ikeuchi 
1988). 

We estimate the mechanical energy input by the starburst 
into the superbubble from the emitted Ha radiation using the 
shock model by Binette et al. (1985). They calculated the ra- 
diative cooling mechanism of shock-heated gas, emitting op- 
tical line radiation, and found that L^ a w 10 _2 L moc i 1 . With 
L Ha = 6.6x 10 40 erg s" 1 for NGC 4410 (MB93) this leads to 



imoch ~ 6x 10 42 erg s _1 . Applying 10 51 erg for the energy re- 
lease per SNll and taking into account that roughly only 20% 
is converted into mechanical luminosity we derive a SNll rate 
of~1.0yr" 1 . 

Under the simplified consideration of a spherically expand- 
ing gas one can estimate its density. With the scaling factor of 
the RS model and setting the electron density equal to the hy- 
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drogen density we get the following expression for the electron 
density of the hot gas: 



2 3D 2 N RS u _ 
nz = — — x 10 cm 



(2) 



where Nrs is the scaling factor, D is the distance to the source, 
r is the radius of the superbubble and / is a filling factor taking 
into account that the hot gas is not distributed homogeneously 
but broken up into separate bubbles. With the parameters of D 
= 97 Mpc, r = 4 kpc, N RS = 5 x 10~ 6 cm 5 and an assumed fill- 
ing factor of 0. 1 one obtains an electron density of 0.03 cm~ 3 . 
Changing the filling factor to / = 0.9 leads to n e = 0.01 cm~ 3 . 

An expansion time of 10 7 yr with a halo gas density of 0.01 
cm~ 3 leads to an X-ray luminosity of the shell of 2.8 x 10 41 erg 
s _1 . A lower SNll rate of 0.5 yr -1 and a slightly smaller ex- 
pansion time of 8 x 10 6 yr reduce the obtained X-ray luminosity 
to 1.2 x 10 41 erg s" 1 in the ROSAT band. The derived plasma 
temperature of 10 7 K lies at the upper bound of the range with 
log T = 6.0 - 6.9 found by S94 for the (0.1-2.2) keV band. 

For a Salpeter IMF, a SNll activity between 10 and 100 
Mq and a SN rate of 0.5 yr -1 the star formation rate results 
to ^95 Mq yr -1 . Depending on the fraction of the mechanical 
energy release of a SNll this value can increase up to a factor 
of 5. 

4.3. Comparison with other galaxies of similar properties 

Each galaxy and, in particular, mergers, galaxy pairs or SB 
galaxies are unique systems. In order to get an insight on 
whether NGC 4410 and its derived structures are somehow typ- 
ical for close encounters, we compare the derived X-ray lumi- 
nosity in the ROSAT band with other disturbed and isolated 
SB galaxies. The peculiar galaxy NGC 2782 e.g. is thought to 
be a merger of two disk galaxies of unequal mass and has Lx 
= 4xl0 41 erg s" 1 (Schulz et al. 1998). Another galaxy with 
disturbed morphology and comparable X-ray luminosity (Lx 
= 1.4x 10 41 erg s- 1 ) is NGC 1808. In contrast the PSPC data 
of this object do not show any X-ray outflow out of the central 
source into the halo. But one has to mention that NGC 1808 has 
a SFR of only 10 M yr" 1 (Junkes et al. 1995). Relatively iso- 
lated systems without any companion, like e.g. NGC 253 (Fab- 
biano et al. 1992), NGC 2903 and NGC 4569 (Junkes et al. in 
preparation), contain X-ray luminosities of a few 10 40 erg s _1 , 
emphasizing the importance of interaction for star-forming ac- 
tivity. 

PC96 found a significant difference in the Lx/Ln a ratio 
between pure AGN, pure SBs and galaxies with circumnu- 
clear star-forming rings with an active nucleus. The pure ac- 
tive nuclei show \og(Lx/Lu a ) between 0.00 and +1.68, while 
the pure SBs in the sample lie between -1.46 and -0.36. The 
three galaxies with combined X-ray emission from AGN and 
SB have values of -0.26 (NGC 1097), +0.16 (NGC 1068) 
and +0.63 (NGC 7469), indicating a continuous decrease 
from AGN to SB. From this tendency one would expect a 
log(L x /X H a)<0 for the RS+PO model for NGC 4410. Our 



results, however yields +0.77 thereby, approximately the same 
as for a single PO model with log(Lx/Lna) = +0.81. 

The fraction of Lh q from the SB relative to the total Ha 
luminosity amounts to 98%, 80% and 40% for NGC 1097, 
NGC 1068 and NGC 7469, respectively. Comparing only the 
contributions from the SB to the Ha and X-ray luminosity, 
PC96 found log(L x /L HQ ) = -0.99, -0.70 and -0.36 for NGC 
1097, NGC 1068 and NGC 7469, respectively. If we concider 
the fraction of Ha luminosity from the SB in these galaxies, 
and assume that 90% of the total Ha luminosity originates from 
the SB within NGC 4410a, L Hq would result in 5.9 x 10 40 erg 
s _1 and log(Lx/ Lua) = +0.36, which is quite high compared 
to the sample analysed by PC96. 



5. Conclusions 

We observed the interacting pair of galaxies NGC 4410 
with the ROSAT HRI and PSPC. Spectral investigations of 
NGC 4410 suggest that the integral X-ray emission (Lx = 
3.9x 10 41 erg s _1 ) can be decomposed into a thermal compo- 
nent (described by a RS spectrum) and a component from the 
AGN (described by a power-law spectrum). The HRI image re- 
veals an extended X-ray halo related to NGC 4410a with an ex- 
tension of 10" from the nucleus of NGC 4410a to the southeast. 
Combining spatial and spectral informations reveals an X-ray 
luminosity of the halo gas of Lx = 1.3 xlO 41 erg s _1 (1/3 of 
the total X-ray emission). The companion galaxy NGC 4410b 
houses only a very faint central point-like source below the 3a 
level, corresponding to an upper limit of 3.8 x 10 39 erg s _1 for 
the X-ray luminosity. 

As a reasonable model we can assume that the tidal inter- 
action in the pair of galaxies NGC 4410 has two effects on the 
one partner, the face-on pec Sab galaxy NGC 4410a: 

(1) A central monster is either formed due to this interaction or 
has already existed before and is now fed by infalling gas dur- 
ing the merging event, producing AGN signatures. Evidence 
for an existing AGN comes from the ROSAT X-ray spectrum 
supported by the spatially correlated radio emission. 

(2) The merging process has ignited a nuclear SB ejecting an X- 
ray gas into the halo. Due to the tidal forces of NGC 4410b, the 
X-ray halo of NGC 4410a is deformed toward its neighbour- 
ing galaxy, either directly by gravitational forces or indirectly, 
originating from a decentred starburst region as a consequence 
of the merging process. 

This system of merging galaxies represents a highly inter- 
esting object because it is a close interacting pair where the 
effects of tidal forces on nuclear activity and SB can be stud- 
ied in detail with large telescopes. The HST WFPC2 image of 
NGC 4410a suggests a possible relation between the deformed 
X-ray halo and a bright emission region near the nucleus. With 
ROSAT imaging it is not possible to distinguish between the 
nucleus and the bright emission region resolved by HST. Fur- 
ther spectroscopic observations of this region and additional in- 
frared imaging could help answering the question of the exact 
origin of the X-ray emission. 
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